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In this work, we theoretically construct exact mappings of many-particle bosonic systems onto 
quantum rotor models. In particular, we analyze the rotor representation of spinor Bose-Einstein 
condensates. In a previous work [I] it was shown that there is an exact mapping of a spin-one 
condensate of fixed particle number with quadratic Zeeman interaction onto a quantum rotor model. 
Since the rotor model has an unbounded spectrum from above, it has many more eigenstates than 
the original bosonic model. Here we show that for each subset of states with fixed spin F z , the 
physical rotor eigenstates are always those with lowest energy. We classify three distinct physical 
limits of the rotor model: the Rabi, Josephson, and Fock regimes. The last regime corresponds to 
a fragmented condensate and is thus not captured by the Bogoliubov theory. We next consider the 
semiclassical limit of the rotor problem and make connections with the quantum wave functions 
through use of the Husimi distribution function. Finally, we describe how to extend the analysis to 
higher-spin systems and derive a rotor model for the spin-two condensate. Theoretical details of the 
rotor mapping are also provided here. 



I. INTRODUCTION 

The behavior of macroscopic systems of multicom- 
ponent bosons under suitable constraints can often be 
greatly simplified through a quantum rotor descrip- 
tion. Within the context of condensed matter physics, 
the most widely appreciated example is the celebrated 
Josephson model [5]. This model provides an accu- 
rate low-energy treatment of two superconductors linked 
by an insulating barrier j2]. The treatment of the full 
many-particle system reduces to a model with two canon- 
ically conjugate variables: the relative particle number 
and phase between the two superconducting regions. 

Bose-Einstein Condensates composed of atoms with 
internal spin or pseudospin degrees of freedom, the so- 
called spinor condensates, offer another arena where such 
rotor mappings are highly useful. Roughly speaking, re- 
cent experimental work investigating the dynamics of 
spinor condensates can be divided into two categories. 
The first category of experiments focuses on the com- 
plex interplay between spatial and spin degrees of free- 
dom resulting from spinor condensates in larger traps 
[SHTO]. These experiments investigate the dynamics af- 
ter a quantum quench, which involves the proliferation of 
topological defects. The second catergory of experiments 
are performed in tight traps where the spatial degrees of 
freedom are unimportant |11H19) . Such experiments have 
focused on the coherent spin dynamics after preparing 
the system in a particular manner. The rotor descrip- 
tion is useful when the spatial degrees of the condensate 
can be neglected, and thus is particularly relevant to the 
second class of experiments. 

In an early theoretical work on spinor condensates it 
was shown that the ground state of the antiferromag- 
netic condensate in tight traps involves large spin cor- 
relations and can be considered to be a condensate of 
singlet pairs of spin-one atoms [20]. However, such "frag- 
mented" states [5TJ are known to be extremely del- 



icate and for most experimental situations are typically 
better described by a broken symmetry state which is 
captured by the classical Gross-Pitaevskii theory [2l5H2"5] . 
Nevertheless the intriguing properties of the fragmented 
condensates in the single-mode regime have motivated a 
considerable amount of further theoretical work 



In this paper we will revisit this problem by employ- 
ing an exact rotor mapping. The mapping, which was 
carried out by some of us in a previous work p], maps 
an antiferromagnetic spin-one condensate in an external 
field onto a quantum rotor model of a particle under an 
external field constrained to the unit sphere PQ. Since 
this mapping is exact, and not a low-energy theory, it 
treats all possible phases of the spin one condensate on 
an equal footing. Roughly speaking, states described by 
the Gross-Pitaevskii Equation (GPE) correspond to ro- 
tor states that are well localized in position. On the 
other hand, states that are delocalized over the sphere 
(e.g. the condensate of singlet pairs of atoms) cannot be 
described by the GPE but are contained within this ro- 
tor treatment. We will provide in-depth analysis of the 
model, and discuss its distinct physical regimes. We will 
also describe its semiclassical limit which has a clearer 
physical interpretation than the GPE and elucidate the 
semiclassical behavior of the rotor wave functions for ap- 
propriate parameter regimes. We will also describe how 
to extend the mapping to systems with larger spin. In 
that sense, the current work is a generalization and ex- 
tension of Ref . [T] . 

The paper is organized as follows. In Sec. |Hj for com- 
pleteness, we consider the simplest nontrivial example of 
bosons in a double-well potential and map this system 
onto a quantum rotor model. We arrive at a result first 
obtained in Ref. [35] but we use a method that can be 
generalized to systems with more components i.e. higher 



spins. In Sec. Ill we move on to the more complex case 
of a spin-one condensate in the single mode regime and 
overview the rotor mapping originally derived by some of 
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us in Ref. pQ. In Sec. IV we consider the correspondence 



between the eigenvalues of the original bosonic problem 
(which has a finite spectrum for fixed particle number) 
and the rotor model (which has an unbounded spectrum 
from above). In Sec. [V] we consider in more detail the 
spectrum of the rotor model, and establish three distinct 
physical limiting cases, namely the Rabi, Josephson, and 
Fock regimes. We provide analytic expressions for the 
low-lying spectrum for these cases 
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we con- 
Here we 



In Sec 

sider the semiclassical limit of the rotor model 
discuss recent experiments on 23 Na dynamics in terms 
the semiclassical phase space. We then connect the quan- 
tum mechanical wave functions in the Rabi and Joseph- 
son regimes to the semiclassical phase space using a gen- 
eralization of the Husimi distribution function 1331 1341 . 



In Sec. VII we consider extending the rotor mapping to 
larger component systems, focusing on the example of the 
spin-two condensates. Finally, in Sec. |VIII| we conclude 
with a summary. 



II. BOSONS IN A DOUBLE- WELL POTENTIAL 

In this Section, we consider the simplest nontrivial case 
of bosons in a double- well potential which is described by 
the dimer Bose-Hubbard model. This archetypal model 
has been studied extensively [SH [351 |3"5H3"T] and has also 
been used to experimentally oberserve the so-called self- 
trapping effect 38J. In the interesting work of Anglin 
et al. [32 it was shown that the dimer Bose-Hubbard 
model can be exactly mapped onto a two-dimensional 
quantum rotor model. Below we will derive their main 
result, using a different formalism which allows a more 
direct generalization to higher dimensional rotor systems 
which will be considered in following sections. We clarify 
our notations and lay out the main theoretical framework 
in this section by considering the double- well case first. 

Our starting point is the Bose-Hubbard dimer model 
which describes bosons in a double-well potential with 
repulsive interactions 



// : = -J(a\a 2 +a\a x ) + ^Un^ni - 1) + ^Un 2 (n 2 - 1). 

(1) 



Here a\ and a\ create bosons in the left and right wells 
respectively, n a = o4a a is the particle number operator, 
J is the hopping, and U is the on-site repulsion. It is often 
instructive to use the amplitude-phase representation of 
the bosonic operators. That is, we can write a — y / ri^e l6 ' Q 
and impose the commutation relation [n a ,#g] = i5 a p. 
Inserting these relations into Eq. (JlJ and expanding to 
leading order in the total particle number N — n\ + 
n 2 (which is taken to be fixed) leads to the well-known 
Josephson model [H |2| 



#jos = -JN cos{d) + Un 2 



(2) 



where 8 = Oi — 8 2 and n — (ni — n 2 )/2 so that the 
two operators in this equation are canonically conjugate: 



[n, 6] — i. The spectrum of the Josephson model can be 
seen to agree with the original double- well model Eq. 
in the large-particle number limit. 

In the work of Anglin et al. |32j . it was shown that 
such a mapping can be made exact, and thus will re- 
produce the spectrum of Eq. ([!]) for arbitrary particle 
number. Their derivation used a method akin to the 
Bargmann phase-space representation of bosonic opera- 
tors [33] ■ Here we will derive their central result through 
a different method. To start, we define the states 
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where fl = (cos(6*), sin(#)) is a real two-component vec- 
tor on the unit circle and the "Cartesian" bosonic op- 
erators b x and b y are defined to be b x = -^{ai + a 2 ), 

by = -^{o-i — 0,2)- These states can be shown to form an 
overcomplete basis. For instance, the fragmented state 
(a|) JV / 2 (a2) Ar / 2 |0) can be seen to be an equal weight 
superposition of these states over the unit circle |22) . 
Therefore, an arbitrary state 1^) in the bosonic Hilbert 
space can be expressed in terms of a superposition over 
the states |Ojv) with weight factor 



I*) = J dn\n N )i;{n) 



(•5) 



Note that due to the overcompleteness, this relation does 
not uniquely determine ip(£l). The approach of the map- 
ping is to find a Hamiltonian H. acting in the "rotor" 
space such that 



dn{H\n N ))i>(n) = / dQp N } (Mip(n)). (6) 



Then the rotor Schrodinger equation Hip = id t ip is a 
sufficient condition for the bosonic Schrodinger equation 
to be satisfied (we will work in units where h ~ 1 unless 
otherwise stated). 

In obtaining "H, we use the gradient operator V on the 
unit circle, which in terms of 9 is V x = — sm(9)dg and 
V a = cos(0)dg. These derivatives satisfy the geometri- 
cally intuitive relations 



V Q 0^ = S, 



a/3 



(7) 



(for a discussion see Appendix [A]) . With this, it can be 
seen that quadratic operators acting on |S~ijv) can be writ- 
ten as 



b%\n N ) = n p (y a + NQ a )\n 



N) 



(8) 



In terms of the Cartesian operators, the double-well 
Hamiltonian (up to a constant offset) is 



H 



U 



j(btb X -blb y ) + ~(ibtby-iblb x f. (9) 
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We can now use the relation in Eq. ^ to find 

{b% - blb y )\n N ) = (iVcos(2#) - sin(20)3 9 ) \fl N ) 



and 



(iblby 



ib\b x 



in 



(10) 



(11) 



where = 



idg. The effective ro- 
tor Hamiltonian can then be obtained by inserting these 
relations into Eq. ([6]), and integrating by parts. One finds 

n = ^Un 2 - J(N + 2) cos(26») - iJsin(26>)n (12) 

where n = idg. While the operator H has a real spectrum 
it is not Hermitian due to its last term. However, one 
can apply a similarity transform to render % Hermitian. 
Specifically, defining [3"2"] 

^ r = e co 8 (2e)^ We -cos(2e)^ ( 13 ) 
and shifting 9 — > 6/2 to compare with Eq. ([2| one finds 



JST = Un 2 - J(N + 1) cos(0) + — sin 2 (0) 



(14) 



which is the main result. Note that this reduces to 
Eq. (pi) in the large- N limit. The additional terms in 
Eq. (125, however, serve to make the spectrum of the orig- 
inal double-well Hamiltonian Eq. ([I]) exactly match the 
eigenstates of this rotor model. 



III. SPIN-ONE CONDENSATES IN THE 
SINGLE-MODE REGIME 

We now move on to discuss the related, but more com- 
plex, problem of the spinor condensate in the single mode 
regime under a magnetic field. Recently it was shown [1] 
that this system maps onto a quantum rotor model under 
an external magnetic field. Here we will summarize this 
mapping. 

Our starting point is a spin-one condensate in a trap 
that is sufficiently tight such that it is a good approxima- 
tion to take all of the bosons to occupy the same spatial 
mode. Under this approximation, we can write the field 
operators for each spin state as 



ip a (r) = (j){r)a a 



(15) 



where a runs from —1 to 1. The condensate profile sat- 
isfies 



d 3 r\4>(r)\ 2 = N 



(16) 



where N is the number of particles in the system. This 
approximation, commonly referred to as the single mode 
approximation, breaks down when the condensate coher- 
ence length is larger than the condensate size. 



The Hamiltonian for this system reads 

u 5 P 2 t 

H = - 1 a o a o- 



(17) 



In this equation, F = a^Fapap is the total spin operator 
where F a p are spin-one matrices, g is the spin-dependent 
interaction, and q is the quadratic Zeeman shift due to 
an external magnetic field. Taking a uniform conden- 
sate density <p(r) = ^/rio, we can express g in terms of 
microscopic parameters as 
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3m 



(a 2 - a )n 



(18) 



where m is the mass of the constituent atoms, and So 
and a,2 are the scattering lengths. We will focus on the 
case of antiferromagnetic interactions for which g > as 
is the case for 23 Na condensates. 

As was done for the double-well problem in Sec. |Hj 
it is useful to transform the bosonic operators to the 
Cartesian basis, rewriting the operators as b x — — (01 — 
a_ 1 )/v / 2, by = (ai +a_i)/iv / 2, and b z = do- Written in 
terms of these, the spin operator becomes F = — iW x b. 
We next define the overcomplete set of states as 



n N ) = 



1 



(n.bt)>) 



(19) 



which are parametrized by a three-component vector on 
the unit sphere O (note that the analogous states in 
Sec.|H] were parameterized on the unit circle). 

The general mapping proceeds with the general 
method given above in Sec.|H] Namely, one writes a gen- 
eral bosonic wave function as a superposition of states 
in the |l~ijv) basis, and finds an operator H acting in 
the rotor Hilbert space which satisfies Eq. ^ (where the 
integration is generalized to the unit sphere). The full 
derivation is given in Ref. [1] and thus we will only give 
the results here. The rotor Hamiltonian corresponding 
to Eq. (fT7l) is 



n = ^l 2 - q (N + 3)n 2 + q n z v z 



(20) 



where L a is the angular momentum operator and V a are 
the gradient operators on the unit sphere. In the spheri- 
cal coordinate representation, V z = — sm(8)dg. This can 
be brought to the more intuitive Hermitian form by ap- 
plying a similarity transformation. In particular, defining 
,^ = e~ s Ue s where S 



2^cos(2(9), we find 



1 



* = 2I L " + V(9) 
where I — N/g is the moment of inertia, and 



(21) 



V(6) = q[N+- )sin 2 (0) + 



q 2 N 



sin 2 (2(9) (22) 



is the external potential. The spectrum of this Hamil- 
tonian exactly matches that of Eq. (17). As was de- 



scribed in [T] one must retain only the eigenstates of 
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this rotor model which are symmetric under inversion: 
-0(17) = ip(—£l). However, since the operator which 
projects into this subspace of states commutes with the 
rotor Hamiltonian, this imposes no additional conceptual 
or technical difficulty. 

Since <fi does not appear in the potential V in Eq. (21 1, 



one notes that this rotor model has azimuthal symmetry. 
This symmetry essentially reduces the model to a one- 
dimensional system which considerably simplifies com- 
putations. One should note, however, that we did not 
rely on this symmetry in the derivation and it will not 
be present for more general couplings. In Appendix |B|wc 
provide a rotor mapping for a more general coupling. 



IV. CORRESPONDENCE OF THE ROTOR 
AND BOSONIC EIGENVALUES 

There are subtleties that arise due to the fact that 



the rotor model Eq. (21) has an unbounded spectrum 
from above, while the spectrum of the original bosonic 
problem for fixed particle number N if finite. As is clear 
from the mapping, an eigenstate of the rotor model ip 
is a sufficient condition for an eigenstate of the bosonic 
Hamiltonian That is, given tp, one can construct the 
bosonic eigenstate through 



*) = / dn\n N )^{n). 



(23) 



Here for simplicity we are taking ip to be an eigenstates of 
the non-Hcrmitian rotor model % so that we do not need 
to include factors of e 5 . Because the spectrum of the 
bosonic Hamiltonian H is bounded, the only possibility 
is that many of the rotor eigenstates get transformed to 
\^f) = through Eq. (23 1, noting that this trivially satis- 
fies the bosonic Schrodinger equation. Following Ref. [32] 
we will refer to the rotor eigenstates which transform 
to | \&) 5^ as "physical" and those that transform to 
|*) = as "unphysical" . 

We next ask if all of the eigenstates of the bosonic 
spectrum are included in the rotor description. For in- 
stance, the pathological case of all the rotor eigenstates 
mapping to \ty) = is not a priori ruled out. Another 
question that arises regards the ordering of the unphys- 
ical and physical eigenstates. In particular, is there an 
energy cutoff below which all eigenstates are physical and 
above which eigenstates are unphysical? We will show 
that there is an affirmative answer to both of these ques- 
tions. 

For sufficiently small particle number N, the eigenspec- 



trum of the bosonic Hamiltonian Eq. (17) can be numer- 
ically computed and compared to the eigenspectrum of 
the rotor system Eq. (21). The rotor Hamiltonian has 



azimuthal symmetry since the potential V appearing in 
Eq. (22) does not depend on the angle tj>. Therefore, 



Jt? commutes with L z and the eigenspectrum for fixed 
values L z = m can be considered separately. Similarly, 
Fx commutes with the bosonic Hamiltonian H, and we 




FIG. 1. Black dashed lines: the eigenvalues E n of the spinor 
bosonic Hamiltonian Eq. ( |17[ ) for N = 12 particles for fixed 
F z — as a function of quadratic Zeeman field q. Red 
solid lines: the lowest 12 eigenvalues of the rotor Hamiltonian 
Eq. (21 1 for fixed L z — after the antisymmetric states for 



which tp(fl) — —ip(Sl) are projected out. The lowest (physi- 
cal) eigenvalues of the rotor model exactly agree with those of 
the bosonic Hamiltonian. For all points in the plot, the spec- 
trum is shifted so that the lowest eigenvalue has zero energy. 



can compare to the rotor model by fixing F z — —m. fn 
Fig. [f] the eigenspectrum of both the rotor model and the 
bosonic model are shown as a function of the quadratic 
Zeeman field q. We take the case of relatively small parti- 
cle number N = 20, and take fixed F z = L z = 0. As can 
be seen, all of the eigenenergies of the bosonic Hamilto- 
nian are accounted for by the rotor model. Furthermore, 
the physical states of the rotor model are always lower in 
energy than the unphysical states. Similar behavior was 
found for other values of fixed L z = m which are not 
shown. 

This behavior can be understood as follows. For un- 



physical states ip{fl), it can be seen from Eq. (23) that 
(Ye m \ip) — for all £ < N. It can be verified that the 
(non-Hermitian) rotor Hamiltonian defined in Eq. ( 20 ) 
has the property 



(Y em \n\Y t , r 







(24) 



for £ < N and £' > N. Suppose that we have a ro- 
tor eigenstate which is unphysical for parameters (q,g). 
Then the eigenstate at (q + Aq, g) can be determined by 
first order perturbation theory. By doing so, one sees 
from Eq. ( 24 ) that if a state is unphysical at q then the 
same state will also be unphysical at q + Aq. In the limit 
of q = 0, the rotor model becomes trivial. Here the eigen- 
states are simply spherical harmonics with eigenenergies 
E( = ^j£{£ + 1). Furthermore, the lowest eigenstates 
for £ < N are all physical while the higher eigenstates 
for £ > N are unphysical in this limit. We note that 
for fixed L z = m, the rotor Hamiltonian becomes one 
dimensional. Thus there will not be any band crossings 
|40) . From the perturbative argument above, we there- 
fore conclude that the higher energy states will always 
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remain unphysical and not mix with the lower energy 
physical states. 



V. EIGENSPECTRA OF THE ROTOR MODEL 

In this Section we will concentrate on the eigenspec- 
trum of the spin-one rotor hamiltonian given in Eq. (21 1. 



We will give the spectrum in particular limiting cases, 
and compare the results with those the Bose-Hubbard 
Dimer problem. 

We consider how the spectrum evolves as a function 
of q. For large q, the potential energy V(0) in Eq. (22) 



serves to localize the wave function on the unit sphere. To 
obtain the energy levels, the angular momentum L 2 can 
be expanded about the north pole so that jtif becomes 
a two-dimensional harmonic oscillator. When 1 <C q/g, 
the second term in the potential energy V{6) dominates 
so that the energy levels are given by 



E (n x ,n y ) = q(n x + n y ) 



(25) 



where n x and n y are integers corresponding to the oscil- 
lator modes in the x and y directions. These eigenstates 
can in fact be directly obtained from the original bosonic 
hamiltonian Eq. (17) in the large-g limit. 



Next we consider the case of smaller q where 1 /N -C 
q/g <C 1. For this case, the first term in the potential 
energy is the most significant. Here we can also expand 
the kinetic energy about the north pole to obtain a har- 
monic oscillator hamiltonian. For this the energy levels 
read 



E, 



/ 2gq(n x + n y ) 



(26) 



of Eq. (17 



where, as in Eq. (25 1, n x and n y are integers. As shown in 



Appendix C]it can be seen that the Bogoliubov spectrum 



agrees with Eqns. (25 26) 



Finally we consider the case of vanishingly small mag- 
netic field such that q/g <§; -^j. For this case the eigen- 
functions are not localized about the north pole. The 
kinetic energy ^jL 2 dominates the rotor model and the 
eigenstates are given simply by 



(27) 



Each of these energy levels has multiplicity 2£+l. So that 
the wave function has inversion symmetry, only even val- 
ues of I should be kept. The ground state in this regime 
where the rotor is completely delocalized about the unit 
sphere, in terms of the bosonic model, is the fragmented 
condensate composed of singlet pairs of bosons. However, 
due to the condition q/g -C 1/N 2 , in the thermodynamic 
limit any small magnetic field will drive the system to a 
symmetry broken state which is described well by mean 
field theory [23 25 . This is the central difficulty in ex- 
perimentally realizing the singlet condensate. We will 
address this problem in more detail in Appendix. [Dj 



It is instructive to compare the above results with the 
dimer Bose-Hubbard model. This model has been an- 
alyzed and found to have three distinct limits, namely 
the "Rabi" , "Josephson" , and "Fock" regimes using the 
terminology of Leggett [23 BTj . Using a method very 
similar to that used above, the expressions for the en- 
ergy eigenstates can be obtained in these regimes from 
the rotor Hamiltonian in Eq. ( Jl4| . Namely, for the Rabi 
regime where N <C J/U the last term in the potential 
energy dominates and the spectrum is approximated by 
a harmonic oscillator, after expanding about 9 — 0. The 
Josephson regime occurs when the first term in the poten- 
tial energy dominates 1/N <C J/U -C N. Here the states 
are also localized about 6 = 0. Finally, for J/U -C 1/N 
the Fock regime is obtained where the rotor is delocalized 
over the unit circle. The Josephson Hamiltonian Eq. ^ 
correctly describes the Fock and Josephson regimes, but 
cannot describe the Rabi regime since a large- ./V expan- 
sion is used to derive it. In summary, the three possible 
regimes for the dimer Bose-Hubbard model are 



A< J/U 
1/N <C J/U -C N 
J/U < 1/N 



Rabi 

Josephson 
Fock. 



(28) 
(29) 
(30) 



It is clear that there is a close parallel between the 
above described regimes for the dimer Bose-Hubbard 
model and those of the spin-one condensate problem. For 
this reason we will adopt the terminology introduced in 
[37l [41] for the spinor problem. Namely, we will label the 
three regimes as 



1 < q/g 

1/N 2 « q/g < 1 
q/g « 1/N 2 



Rabi 

Josephson 
Fock. 



(31) 
(32) 

(33) 



For typical experimental situations (for example those 
describe in Refs. [HI Qjj]) 9 ~ 9 and N ~ 10 3 - 10 5 
which places the system in either the Rabi or Josephson 
regimes. For these cases, the Gross-Pitaevskii equation 
gives a qualitatively correct description of the dynamics. 
It is also interesting to note that for double-well conden- 
sates the Rabi regime is more difficult to achieve since 
by reducing the hopping J to achieve the condition in 
Eq. (28), a single-band description becomes inapplicable. 



On the other hand, the Fock regime for double-well con- 
densates can be experimentally achieved, which has the 
Mott Insulating ground state [HJES]. 



VI. SEMICLASSICAL ANALYSIS OF THE 
ROTOR MODEL 

In this Section we analyze the semiclassical limit of 
the rotor model in Eq. (21). We will use this to address 
recent experimental results. We will then show results 
from taking the Husimi transform of the quantum eigen- 
states of the rotor model. Such a method has been shown 
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FIG. 2. Left: Equal-energy contours of the semiclassical energy given in Eq. (351 for 



= 0. Right: Husimi distribution 
functions H(z) for particular eigenstates of the rotor model Eq. (21 1 for the parameters q = g, N = 10, and « = 1/10. Panels 
(b), (c), (d), (e) respectively correspond to the 2nd, 4th, 6th, and 8th excited states within the manifold of m — and even 
I. The classical equal-energy contours corresponding to the energies of the states plotted on the right are shown in red. The 
same range of pg and is used for all plots. 



to elucidate the semiclassical limit of the Bose-Hubbard 
dimer model [34] , 

In recent experiments [HI [T!5] the dynamics of a 23 Na 
condensate, which has antiferromagnetic interactions g > 
0, was investigated. The initial condensate was prepared 
in a fully polarized ferromagnetic state pointing in the x 
direction after which the condensate was allowed to freely 
evolve. The value of (F x ) was measured as a function 
of time. Two distinct types of behavior were found, de- 
pending on the external magnetic field which couples to 
the system through the quadratic Zeeman shift q. A sep- 
aratrix between these two behaviors occurs at a critical 
magnetic field B c . When B < B c , (Fx) 2 showed oscil- 
latory behavior, having (Fx) 2 > at all times. On the 
other hand, when B > B c it was seen that (F x ) 2 = 
at periodic intervals during its evolution. An analysis of 
the behavior was provided in terms of the classical Gross- 
Pitaevskii energy functional. Taking into account the 
conserved quantities (total particle number and spin mo- 
ment in the z-direction which was fixed to be (F z ) = 0), 
the two-dimensional phase portrait of the energy func- 
tional was shown to capture these two regimes. 

We will now describe how the semiclassical limit of 



Eq. (21 1 gives an intuitive understanding of these results. 



The corresponding Lagrangian is 



C 



+ sin 2 (6»)0 2 -V(0) 



where V(0) is given by Eq. (22 1. This gives the canonical 
momenta pg = 19 and p^ = Isxn (ff)<j>. The correspond- 
ing classical energy is 



E = 




(35) 



The above equations describe the motion of a particle 
on a unit sphere. We will concentrate on the case where 
P0 = 0. The classical equal-energy contours of Eq. (35) 
are plotted in the left panel of Fig. [2] Two types of 
behavior are seen. The higher-energy states have motion 
where the particle's trajectory explores both hemispheres 
but has either pg > or pg < at all times thus never 
having zero angular momentum. This corresponds to the 
motion of the spin-one condensate for B < B c . Increasing 
the magnetic field will constrain the particle's trajectory 
to one hemisphere. For this motion, it is seen that pg = 
at periodic intervals during the particle's trajectory. This 
type of motion corresponds to B > B c of the spin-one 
condensate experiment [4"4"| . 

We will now move on to discuss the semiclassical prop- 
erties manifest in the quantum mechanical wave func- 
tions of Eq. (21) for appropriate parameter regimes. To 



(34) 



do this, we will use a generalization of the Husimi dis- 
tribution function (33] to the case of the sphere. The 
Husimi distribution function has been successfully ap- 
plied to elucidate the pendulum structure manifest in the 
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wave functions of double-well condensates [51] described 
by Eq. 0. 

Conventionally, the Husimi distribution function of a 
particular wave function \ip) is defined as 



H(z) 



W>| 2 



(36) 



where |z) is a scaled coherent state. For our consider- 
ations, we thus need a generalization of the notion of a 
coherent state to the unit sphere. Recent work on such a 
generalization is given in Refs. [13J HB] . In [IS] it was ar- 
gued that it is most natural to define spherical coherent 
states to be eigenstates of the "annihilation" operators 



n, 



1 t- T 2 
P 2 KL 



(37) 



where k is a scaling parameter (47) . Such eigenstates are 
given by 



\ Z )=J2e~ Kt{e+1)/2 \Yt m )Y; m (z) 



(38) 



I'm 



In this equation z is a three-component complex vector 
that satisfies z • z = 1. In terms of classical phase space 
variables (p = p e 9 + p,p(f> and O), z can be expressed as 



z = cosh(ftp)f2 +i— sinh(/tp)p. (39) 
P 

The value of the scaling parameter should be taken such 
that k 2 ~ ^2 which is the ratio of the prefactors of the 



kinetic and potential energies in Eq. (21 1. 

We consider the case where q — g and N = 10 bosons 
which places the system between the Josephson and Rabi 
regimes, and well away from the Fock regime. Density 
plots of the Husimi distribution function for particular 
rotor eigenstates are shown in the right of Fig. [2] Since 
we concentrate on the case of F z = 0, the Husimi distri- 
bution function will only depend on the pair of variables 
(6,po). Red classical equal-energy contours shown on 
the right of Fig. [2] are drawn for energies corresponding 
to these eigenstates. One sees that the Husimi distri- 
bution functions strongly resemble the semiclassical con- 
tours, thus revealing the semiclassical behavior of the 
eigenstates. Such agreement occurs for all states in the 
Rabi and Josephson regimes, but not for the Fock regime 
which has no semiclassical correspondence. 



VII. EXTENSION TO HIGHER SPINS 

We will now move on to discuss how to perform the 
rotor mapping for larger spin systems. We will focus on 
F = 2 spinor condensates because of their experimental 
availability as hyperfine states of alkali atoms. We will 
show that this system maps to a particle moving on a 
sphere in five-dimensional space. 



Spin-two condensates have five spin components. We 
take a a for a — —2 ... 2 to annihilate, a boson with 
F z = a. In the single-mode regime, spin-two condensates 
are described by the Hamiltonian [JS] 



H=^F 2 
2N 



2N 



(40) 



Here, F = a^ a F a pap is the total spin operator where F a p 
are spin-two matrices. In the second term, A = aoao — 
2a 1 a_ 1 + 2a 1 a,-\ annihilates a singlet pair of bosons. In 
terms of physical quantities, the coefficients g\ 2 are given 
by 



9l 



.92 



4tt?i 
7m 



n (a 4 - 02) 



n 



5 («o 



a a) 



-(o 2 - a 4 ) 



(41) 
(42) 



where ao, a.2, and S4 are the spin- two s-wave scattering 
lengths. For simplicity we will neglect the effects of an 
external magnetic field on this system. 

We perform the following unitary transformation on 
the bosonic operators: 



h 
hi 

h 

b 4 

h 



ao 
1 

iy/2 
1 

1 

iV2 
1 

T2 



(-a_i - ai) 
(a_i - 01) 
(02 - a_ 2 ) 
(a 2 + a_ 2 )- 



(43) 
(44) 

(45) 

(46) 

(47) 



These operators transform as a vector under SO (5) rota- 
tions generated by M a p = —i(b\j)p — bpb a ). In terms of 
these quantities, the singlet operator is 



^ = bb 

while the spin operators are 

F x = V3M 12 - M 25 
F y = V3M 13 + M 24 
F z = M 23 + 2M 45 . 



M 34 
M35 



(48) 



(49) 

(50) 
(51) 



As before, we can parametrize an overcomplete set of 
states (but now using the five- component, real unit vector 
fi) as 




(52) 



The mapping proceeds along similar lines to that in 
Sees. [Tl] and |III| In particular, one finds that 



M, 



a/3 



-L 



a0 



(53) 



where L 



a/3 



-i(Q a V£ 



fl^V a ). This can be used to 



find the rotor correspondence of the first term in Eq. ( 40 1 



Next we find the rotor correspondence of the second 
term in Eq. ( 40 1 . We use the five-component version of 
Eq. (|8| to find that 



A*A\fl N ) = (V 2 + N 2 + 3JV) |n 



(54) 



In this equation, V 2 is the Laplacian on the five- 
dimensional hypersphere, as described in Appendix [A] 
The integration-by-parts here is straightforward. One 
obtains 



A* A -> V 2 + N 2 + 3N. 



(55) 



The resulting rotor model, which is already Hermitian, 
is thus 



* =^ V 2 + |L ( (V3L 12 - L 25 + L 34 ) 2 + ( r, (i ) 



2N 

(\/3L 13 + L 2i + L 35 ) 2 + {L 23 + 2L 45 ) : 



where we have dropped a constant energy offset. This 
model has a particularly simple form in the limit of g\ = 
0. Here the system has an SO (5) symmetry, and the 
ground state will be a condensate of singlet pairs of spin- 
two bosons. 



VIII. CONCLUSION 

In this work we have analyzed in detail rotor mappings 
of spinor condensates in the single mode regime. We have 
addressed some subtleties related to the physical and un- 
physical eigenstates and showed that the rotor mapping 
gives an exact treatment of the spinor condensate. Since 
the rotor model treats the mean field as well as corre- 
lated phases on equal footing it offers new insights into 
the problem. We have established both the importance 
of the rotor model in providing physical insight into the 
properties of spinor condensates and its validity as a prac- 
tical scheme for carrying out calculations. 

There are several interesting directions that can be 
pursued in future work. The Husimi distribution func- 
tion has proven useful for understanding the collapse and 
revival process of atoms in the proximity of a superfluid- 
insulating phase transition [43] . Such a phase-space anal- 
ysis of the collapse and revival dynamics for the spinor 
system close to the Fock regime will prove to be valuable. 
We emphasize that for this regime the semiclassical cor- 
respondence illustrated in Fig. [2] will not hold. 

In Sec. |VII] we derived the rotor representation of the 
spin-two system for a single site. The mean-field phase 
diagram of the spin-two condensates is known to have 
a degeneracy for nematic states [JS] which is lifted by 
quantum and thermal fluctuations [SOI HI] • A general- 
ization Eq. ( 56 1 to include quadratic Zeeman field will 



prove useful for studying this effect for smaller conden- 
sates where quantum effects are more pronounced. Fi- 
nally, we note that low energy effective rotor theories of 



spinor condensates (without magnetic fields) where pre- 
viously investigated in 52-54 . It will be interesting to 
investigate how the rotor mapping generalizes to include 
spatial degrees of freedom. 



ACKNOWLEDGMENTS 

It is our pleasure to acknowledge useful discussions 
with A. Lamacraft, P. Lett, and W. Reinhardt. This 
work was supported by the NSF Joint Quantum Insti- 
tute Physics Frontier Center. We thank the hospitality 
of the Kavli Institute of Theoretical Physics under the 
grant NSF PHY05-51164 where part of this work was 
completed. 



Appendix A: Quantum mechanics on the 
hypersphere 

In this Appendix, for convenience, we will tabulate the 
properties of a particle constrained to the surface of a d- 
dimensional hypersphere. The position of the particle 
is given by d coordinates fl a (for a = l,...,d) sub- 
ject to the constraint O • f2 = 1. The momentum op- 
erators are ir a = — iV a where V a is the gradient oper- 
ator in the a-direction on the hypersphere (which can 
be expressed in terms of d — 1 angles and their deriva- 
tives). Finally, the angular momentum operators are 
L a p = fl a Ttp — QpTTa. Note that for d = 3, the an- 
gular momentum is conventionally written with a single 
subscript as L a — ie a p^Lp^. The position and angular 
momentum are Hermitian operators, while the Hermitian 
conjugate of 7r Q is 



tt„ = 7T Q + i(d - 1)0 Q 



(Al) 



The following commutation relations are satisfied for the 
position and momentum operators: 



[n a ,n p ] = o 

[Q a ,Trp] = i(5 a p 



-iL, 



a/3- 



(A2) 
(A3) 
(A4) 



These can be seen to give the angular momentum op- 
erators the correct commutation relations which are 

[L a p,L~/s] — iSayLps + iSpsLay — iSagLp^ — iS/3jL a S- 

These operators satisfy the orthogonality relation O • ir = 
(i.e. the momentum and position are orthogonal on the 
hypersphere) . It can also be verified that the total angu- 
lar momentum can be expressed as 



(A5) 



Appendix B: Rotor model with general coupling 

In this Appendix we consider spin-one Hamiltonians 
with more general coupling to external fields. In partic- 
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ular, we consider Hamiltonians of the form 



H 



2N 



,F 2 + biB a0 bp. 



(Bl) 



where B is a Hermitian matrix. One can see that this 
reduces to Eq. (17) for the special case B a p = —q5 az 6p z . 



In the following it is useful to write B in terms of its 
real and imaginary parts as B = B' + iB" . Since B is 
Hermitian, we have the requirement that B' is symmetric 
while B" is antisymmetric. 

Applying the rotor mapping as in Sec. |III| one arrives 
at the non-Hermitian Hamiltonian 



It is straightforward to diagonalize this by a Bogoliubov 
transformation. The result is 



H = y/q(2g + q)(^a + /3 f /3) (C3) 
where a and /3 are bosonic annihilation operators which 
is consistent with |30j . The spectrum of this hamiltonian 
can be seen to agree with the results derived from the 
rotor model in Sec. [V] in the Rabi regime, Eq. (25), and 
Josephson regime, Eq. (p6|. However, the results do not 



agree in the Fock regime, Eq. (27), since the Bogoliubov 



treatment is inapplicable for a fragmented condensate. 



H 



2N 



L 2 + q(N + 3)B aP Q, a n p - B aP n p V a (B2) 



which should be compared to Eq. (20). To bring this 



Hamiltonian to Hermitian form we apply the similarity 
transformation = e~ s T-le s where 



T a pfl a flp 



(B3) 



where T is a matrix. One can verify that by choosing 
r = jjjB' , provided B' and B" commute, Jtff becomes 
Hermitian. In particular, for this value of T, 



J?=^ I L 2 + V(9,0) 



where 



v(e,<t>) = ( n+ | n T B'n 



N 
2^ 



n T B' 2 n - {n T B'ny 



(B4) 

(B5) 
(B6) 



One can check that this reduces to Eq. (21) in the ap- 
propriate limit. 



Appendix C: The Bogoliubov spectrum of Eq. (17) 



It is instructive to compute the low lying spectrum of 



the spinor Hamiltonian Eq. (17) through the Bogoliubov 
method [SS] and compare with the results from the exact 
rotor mapping given in Sec.[Vj We expand about classical 
mean- field state given by a\ = a_i =0 and ao = yN, 
and write the bosonic operators as a a — a a + 5a a . The 
constraint of fixed particle number iV can be enforced up 
to quadratic order by requiring 



-y/N(6. 



a 



5a 



5a[Sai + SaL^a-i. 



(CI) 



Dropping constant terms, Eq. (17) becomes up to 
quadratic order 

H = (g + q)(8a\Sai + 6a_i8a-i) + g{5a\5a^i + h.c). 

(C2) 



Appendix D: Experimental realization of the singlet 
condensate 



In this Appendix, we will discuss the experimental pa- 
rameters necessary to achieve the singlet condensate. In 
the Josephson regime, the ground state wave function of 
the rotor model is 



(Dl) 



where 6 = 



As can be verified from ( 32 ) , in the 



Josephson regime, « 1. Decreasing the magnetic field 
and thereby decreasing q, one sees that the width of the 
wave function increases. When the width of the wave 
function 9 approaches unity, the harmonic description of 



the condensate fails and the Fock regime is entered ( 33 1 . 
As mentioned earlier, in the limiting case of q = the 
ground state is a condensate of singlet pairs of spin one 
bosons. 

We thus ask what parameters are necessary for 9 ~ 1. 
Because this quantity scales inversely with the number 
of particles, this state cannot be achieved in the ther- 
modynamic limit. We therefore concentrate on systems 
with relatively small particle number. For the 23 Na sys- 
tem, the parameters g and q appearing in the rotor model 
Eq. (21) are related to the atomic density and external 
magnetic field B through [18] 



g = (1.59 x l(T 52 Jm 3 ) n a 

q = (1.84 x KT 35 J/(//r) 2 ) B 2 



(D2) 
(D3) 



For a fixed particle number, increasing the density in- 
creases 9. We thus take no = 5 x 10 14 cm -3 , which is 
relatively large, but still small enough that three-body 
losses are not important. Then for small magnetic field 
B = O.l^iT and N = 500 particles, we have 9 = 1.9 which 
is outside of the Josephson regime. If quenched from fi- 
nite magnetic field, such a system will exhibit quantum 
collapse and revival oscillations [T]. 
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